Nonvanishing spin Hall currents in the presence of magnetic impurities 
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The intrinsic spin Hall conductivity in a two dimensional electron gas with Rashba spin-orbit 
coupling is evaluated by taking account of anisotropic coupling between magnetic impurities and 
itinerant electrons. In our calculation Kubo's linear response formalism is employed and the vertex 
correction is considered. In the semiclassical limit /j, ^> 1/r, a non- vanishing spin Hall conductiv- 
ity o sH is found to depend on the momentum relaxation time r, spin-orbit splitting A and the 
anisotropic coefficient of interaction between itinerant electrons and magnetic impurities. The clean 
limit of a aH is in the region of e/8n ~ e/6n, depending on the anisotropic coefficient. 

PACS numbers: 72.25.Dc, 72.25.-b, 72.25.Rb 
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I. INTRODUCTION 

Recently the spin Hall effect (SHE) [ESSSISIEBU 
has attracted much attention due to its potential applica- 
tion in spintronics. In the spin Hall effect, a longitudinal 
electric field creates a transverse motion of spins with the 
spin-up and spin-down carriers moving in opposite direc- 
tions, which leads to a transverse spin current perpen- 
dicular to the external electric field. To understand the 
spin Hall effect, one needs to study the intrinsic spin Hall 
effect (ISHE) that has been discussed intensively. Theo- 
retically ISHE may exist in the p-type semiconductor [l[ 
and two-dimensional electron gas (2DEG) 0]. Sinova 
et al. predicted a universal spin Hall conductivity in 
clean 2DEG as 



° sH = 8^ 

From then on, many researchers have devoted to the issue 
of whether this result can be modified in the presence of 
impurity. Some authors found that an arbitrarily small 
concentration of impurities would suppress the spin Hall 
conductivity to zero due to the vertex corrections 

SOS 

EHH El- While others 0H El argued that the spin 
Hall conductivity was robust in the presence of disorder, 
falling to zero only when the lifetime broadening is larger 
than the spin orbit splitting of the bands, i.e., 1/r > A. 
Grimaldi [17| considered the situation of sufficiently low 
electron density, and found that the vertex corrections 
no longer suppressed the spin Hall conductivity when the 
Fermi energy is comparable to or smaller than the spin- 
orbit energy. Other authors revealed that the vanishing 
of <j s h was a peculiar feature of the linear Rashba model. 
Taking into account of nonlinear momentum dependence 
of the spin-orbit interaction a(p) [hsj or a nonquadratic 
band spectrum e(p) [19[, u s h was robust against impurity 
scattering. 

Very recently, an important step was made by Inoue 
et al. (20| who found the spin Hall conductivity was not 
zero in the existence of magnetic impurities in the limit of 
At ^> 1. Physically, the acceleration of the electrons by 



the external electric field modifies the SO-induced pseu- 
domagnetic field such that the spins are tilted out of the 
2DEG plane in directions that are opposite for positive 
and negative lateral momentum states. This corresponds 
to a flow of a z = 1/2 and <r z = —1/2 spins in oppo- 
site directions without a corresponding net charge trans- 
port. In the presence of isotropic impurity scattering, the 
spin Hall current is proportional to the time derivative of 
the spin polarization [13j | which vanishes in a stationary 
state. Whereas, this relation does not fulfil for magnetic 
impurities, leaving over an opportunity for nonvanishing 
spin Hall conductivity. As we are aware, the case con- 
sidered by Inoue has not been developed to anisotropic 
model in the semiclassical limit fi 3> 1/r. And Inoue 
et al. adopted a simplified approximation in calculating 
the integral of Green's functions. It is therefore obliged 
to study such kind models and investigate spin Hall con- 
ductivity carefully. 

In this paper we calculate the spin Hall conductivity 
for a two-dimensional magnetically disordered Rashba- 
electron gas where the magnetic interaction is anisotropic 
in the limit of large Fermi circle fj, ^> 1/r, A. The XXZ- 
type interaction between the magnetic impurity and the 
electron spin is adopted. Our calculation is carried out by 
considering the vertex corrections with the help of Kubo's 
linear response formalism within the self-consistent Born 
approximation. The paper is organized as follows: In 
Sec. mi we introduce the model by taking into account 
of magnetic impurities. In Sec. IIII1 the calculation pro- 
cedure of the spin Hall conductivity is presented. In 
Sec. lIVl we give a comparison between our result and the 
result in other disordered system and discuss the possible 
reasons that cause the vanishing spin Hall conductivity. 



II. THE MODEL WITH MAGNETIC IMPURITY 

We consider the 2DEG with Rashba spin-orbit cou- 
pling in the presence of impurities, whose Hamiltonian 
consists of two parts H = H + Vdis with H the sum 
of the kinetic and Rashba terms, and Vdis the potential 



caused by impurities. The main contribution of a short- 
ranged magnetic impurity can be described by an inter- 
action between an itinerant electron and a local magnetic 
moment (2l| whose orientation is defined by polar and az- 
imuthal angles (9,4>) (see Fig. Q}. A quite general form 
of such kind interaction is of XXZ type whose second 
quantization form is given by the following Hamiltonian 



N 

Vdis = z2 dr 2 uS{r - R, 



(1) 



7COsWi smff.e 
sin#je l< " — 7Cos( 



where (Ri,#i,</>i) denote the position and orientation of 
the ith impurity N the total number of impurities; u and 
7 refer to, respectively, the strength and the anisotropy 
of the coupling between the itinerant electrons and the 
impurities. 



1>l (r) 



(2) 



The interaction becomes isotropic for the special value 
7 = 1 which corresponds to the model considered by In- 
oue We assume that the distribution of impurities 
is homogeneous and their orientations are isotropic, ac- 
cordingly, 



P(Ri0i<£i,R 2 02</>2,--- ,~R N d N( f> N ) 
= ( — ) N dR 1 dn 1 dR 2 dn 2 ■ ■ ■ dR N dil 



^ nV < ,—, — »--n, (3) 

where dO,{ = sin B ^dB%d(j)i and V denotes the area of 
the 2DEG. Note that the model given by Eq. ((TJ for 
magnetic impurity is not reducible to that for nonmag- 
netic impurity since the matrix between ijj' (r) and ip(r) 
for the former is a traceless matrix while that for the 
latter is an unit matrix. For a general discussion, the 
strength of Rashba coupling is momentum-dependent 
i.e., H = e(p) + a(p)(a x Py — <TyPx). This Hamiltonian 
Hq can be diagonalized to be 



e±(p) = e(p)TJJa(p), 



by the unitary matrix 



U(P) 



1 
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(4) 



(5) 



where ip p denotes the angle between the momentum p 
and the x axis. The momentum-dependent Rashba cou- 
pling has important effects on the spin-Hall conductivity 
in the presence of nonmagnetic scatterers [18| , while it is 
not relevant to the final result in the present case with 
magnetic scatterers. 



III. SPIN HALL CONDUCTIVITY 

As is well known that spin angular momentum is not 
conserved in the presence of spin-orbit coupling, and the 




FIG. 1: Bloch sphere representation of an magnetic impurity. 



spin current is defined as 



1 



{Vy<J z + <T z Vy) 



4 

1 de(p) 

2 dp y 



(6) 



On the basis of Kubo's formalism, the spin Hall conduc- 
tivity can be expressed as [l3| 

asH ^ = ^v J l^ Tr {" F(cj + Wl)J y 

[G r (cu + ui) - G<> + wi)] jxG a {u x ) 

+ n F (u Jl )J^G r (u J + uj 1 )j x [G r (u Jl )-G a (uj 1 )]Y (7) 

where n F (u>) — l/(e /3l - UJ ~^ + 1) is the Fermi function, 

jx{p) = -r~[s(p) + oi(p)(a x py — a v p x )] is the one-particle 
dp x 

charge current operator, G r and G a are the retarded and 
advanced Green's functions. The trace is taken over mo- 
mentum and spin indices. The diagrammatic method is 
employed to calculate the average spin Hall conductivity 
over the distribution of impurities. The trace in Eq ■([?]) 
is expanded into the sum of diagrams (see Fig. 

A. Average Green's functions 

Now we calculate the average Green's functions in the 
self-consistent Born approximation. For simplicity, the 
Fermi energy /i is set to be the zero point of energy since 
it can be combined in the dispersion relation. Then the 
bottom of conduction band goes to negative infinite in 
the large Fermi circle limit. The free Green's function in 
chiral bases can be expressed as 



1 



G 0(ch)(P^) = 



w-e+(p) + iv 




-e-(p) +ir) 



(8) 



where s±(p) — e(p) =p pa(p). The Feynman diagram 
for the magnetic impurities is shown in Fig. [3] The 
self-consistent Born equation for the averaged retarded 
Green's function in chiral bases is given by 



3 



Tr [J*G r j x G a ] = 




G r G r G r 



FIG. 2: The spin Hall conductivity is calculated in terms of diagrams. The first diagram contributes to cr sH , while the others 
contribute to the ladder correction for spin Hall conductivity. G r and G a denote the average Green's functions taking over the 
distribution of impurities. 
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FIG. 3: The Feynman diagram for the magnetic impurity (left hand side of the diagram equation) is different from that for 
the nonmagnetic impurity (right hand side). Each dot line connected to the same impurity gives an extra matrix and the 
integration J d6d<j> sin 9/4-k for 9 and 4> is performed. 



q ch) (p,a;) = GS (c , ) (p )W ) + W d ed^^Gl (ch) ^u)uH^)(Z°e^ -7^)^ 

q 

G( ch) (q,a;)[/t (q) ( 7™* ) [/(p)G[ c/l) (p, u ), (9) 

I 

The Green's functions in a z bases are 



G a (p, uj) = U(p)G a ch (p, w)U\p). (12) 



Then we calculate o° sH , which corresponds to the one- 
loop diagram in Fig. 2, at zero temperature and zero 
frequency 




where N denotes the number of impurities. Eq. ([9| has 
a solution 



1 



G (cfc)(P> W ) 



u-e+(p) + 27 




(10) 



where r denotes the momentum-relaxation time and 
1/t = 7i 1 u 2 2ttN f {-/ 2 + 2)/3, Hi = iV/V the impurity 
concentration and Np the density of states at Fermi sur- 
face. Similarly, the averaged advanced Green's function 
is 



1 



w-e+(p) ~ 27 




w-e_(p) - 27 / 



(11) 



4 



sH 



Um ±- I ^ ^Tr{0(- W - Wl )^(p)[G r (p, w + Wl ) - G a (p, w + wi)]i s (p)G°(p, "0 
p 

+ 0(- Wl ) ^(p^p, w + Wl )i x (p)[G r (p, Wl ) - G a (p, Wi)] } 



1 



1 + (At) 5 



(13) 



r 



where A = 2p F a(pp) is the spin-orbit splitting at the 
Fermi surface, and p F | £ ( p is the Fermi momentum. 
In deriving the last line of Eq. (| 13[> . we took the integral 
with respect to uj± and the summation with respect to p 
in the limit of large Fermi circle /j> 1/t, A by means of 
the same method in Ref. fl3[. C learly, a® H is the same as 
that derived by Dimitrova [131 ] for nonmagnetic impuri- 
ties. But it differs from e/8n derived by Inoue [20|, who 
takes the product of Green functions as a 5-function for 
simplicity. Such an approximation is valid when At S> 1, 
but it is inappropriate to take this approximation before 
taking account of the vertex corrections to the spin Hall 
conductivity. 



correction with one advanced and one retarded Green's 
functions is considered. The sum of ladder diagrams at 
to = gives 



'sH 



2irV 



Y,Tr[J*G r (p,0)j x (p)G a (p,0)], (14) 



where the sum of the series of vertex corrections to the 
Jy(p) is denoted by the matrix J* (see Fig. 0]). 



B. Vertex correction 

Furthermore, we calculate the vertex correction a^ H 
to the spin Hall conductivity, which corresponds to the 
ladder diagram in Fig. 2. Vertex corrections to the terms 
in Eq.([7]) with two advanced or two retarded Green's 
functions vanish as fx ^> 1/t, A [13]. Only the vertex 



+ 




FIG. 4: The vertex of spin current with vertex corrections. 



The self-consistent equation of J* is 



riiU 



E 



1 . „ ( j cos 9 sinOe l ^ 
4^ Isinfle** -7COs6> 



G a (p,0)(J*(p) + J*)G r (p,0) 



7 cost/ smf/e 



sin ( 



— %4> 



-7 cos ( 



(15) 



r 



In the limit of large Fermi circle ijl 3> A, 1/t, the sum- 
mation over momentum in Eq. (|15p can be evaluated by 
taking integral. As a result, we have 



(^)tt 








7 2 

yG)(J*) u , 






yG)(J*) TT 4 






(Jy)n 


7 2 


- y5(J») Ti) 






( J y)ll 


I 2 

=v- 


7 2 ~z 

-J B { J y)l\i 




(16) 



3iv F Ar 



and advanced Green functions 
A = 
B = 
C = 



4( 7 2 + 2)[1 + (At) 2 ] 
1 



3 .1 



7 2 + 2 v 2 2[1 + (At) 2 ] 
1 



3 -(- 



7 2 + 2 v 2 2[1 + (At) 2 ] 
From these equations, we find (Jy)fi = (Jy)u and 



(17) 



yJU 



<Jy) 



yM 



iv F AT 



2+(8/ 7 2 + 6)[1 + (At) 2 



(18) 



where v F = de(p)/dp \ e u,\ =li is the Fermi velocity. A, B, 

and C are momentum integrals over products of retarded The vertex correction to the spin Hall conductivity is 
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evaluated as 



'sH 



-ieAr(J') 



y'U 



1 



4:7TVp 1 + (At) 2 

(At) 2 



l + (4/ 7 2 + 3)[l + (Ar) 



(19) 



Summing Eq. (fTU)) and Eq. lfT!))) . we obtain the spin Hall 
conductivity 



<?sH 



(At) 



'sH 



J sH 



8tt 3 7 2 
1 + 



(20) 



4 7 2 



"(At) 2 



This is our main conclusion. 

Now we discuss the relation between our work and a 
recent work on such topics. Inoue et al. [20J considered 
the vertex corrections due to magnetic impurities, which 
correspond to the isotropic case 7 = 1 of our model. In 
this special case, our result becomes 



CsH = (1 



8tt 



(At) 2 



(1 



f 7(At) 2 ' 
(At) 2 
8 + 7(At) : 



sH 



r)d- r 



(At): 



:)• (21) 



This differs from the result obtained by Inoue et al. 2 
In their paper, the product of Green functions in Eq.fTSl 



was taken to be a (^-function so that their a® H and a^ H 



lack the factor (1 



1 



Such an approxima- 



l + (Ar) 2 ^ 

tion does not affect the total spin Hall conductivity in 
the existence of nonmagnetic impurities because the a aH 
and <Tg H cancel each other. However, the cr° sH and <j^ h 
do not cancel each other in the magnetically disordered 

system, and hence the factor (1 — — -tt) enters the 

J v 1 + (At) 2 ' 

final result. Even in the situation of tA ^> 1, the ap- 
proximation condition of Ref. [2(1, our ct s //| tA>1 docs 
not have the same lower order approximation as that 
e „ (At) 2 



8^ (1 "' 8 + 7(At) : 



-) in Ref. |20|. 



a Rashba-split 2DEG. It was assumed that the magnetic 
impurities were short-ranged and the orientations of their 
local moments were distributed uniformly. Unlike the 
case of nonmagnetic impurities, the vertex correction a^ H 
does not cancel <J® H , leading to a non- vanishing spin Hall 
conductivity (see Eq. (f20"]) ) which depends upon the spin- 
orbit splitting A and the momentum relaxation time r in 
the limit of large Fermi circle /j > A,l /t. The scattering 
changes both the momentum and spin directions from a 
magnetic impurity and lead to some correlation between 
them. Therefore, the average spin tilting cannot cancel 
the spin current completely. In the dirty limit t <C 1/A, 
&sH goes to zero. While in the clean limit r — ► 00, the 
spin Hall conductivity is 



lim a sH 



e 4 7 2 + 4 
8^3 7 2 + 4' 



(22) 



which depends on the coefficient 7. When 7 changes 
from to 00, the clean limit of er s # changes from e/87r 
to e/67r, and is e/77r for the isotropic interaction. This 
result differs from the universal value e/87r for the ideal 
clean 2DEG and also differs from the result a s n = in 
nonmagnetically disordered system. 

Concerning to the vertex corrections, the spin Hall 
conductivity in nonmagnetically disordered or magnet- 
ically disordered systems does not go back to the uni- 
versal value e/87r even if the impurity concentration is 
arbitrarily small. This result comes from the infinite size 
of system and infinite decoherence length L v assumed in 
calculation. As r — > 00, the corrected vertex (Jy)-\i goes 
to zero, but the summation of Green's functions diverges 
because of lack of a cut-off in momentum space. This 
divergence leads to a non- vanishing a^ H and changes the 
total spin Hall conductivity. In real systems, the quan- 
tum interference contributing to the vertex corrections 
only happens inside the decoherence length. So the spin 
Hall conductivity depends on the decoherence time t^ 
when the momentum-relaxation time is large enough to 
be comparable with it r ~ t v . The a S H derived in this 
work is only valid when r is not too large. Further studies 
are in progress. 



IV. SUMMARY AND DISCUSSION 
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